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ABSTRACT 

We present dynamical distance estimates for 15 Galactic globular clusters and use these to check the con¬ 
sistency of dynamical and photometric distance estimates. For most of the clusters, this is the first dynamical 
distance estimate ever determined. We extract proper-motion dispersion profiles using cleaned samples of 
bright stars from the Hubble Space Telescope proper-motion catalogs recently presented in Bellini et al. (2014) 
and compile a set of line-of-sight velocity-dispersion profiles from a variety of literature sources. Distances 
are then estimated by fitting spherical, non-rotating, isotropic, constant mass-to-light (M/L) dynamical mod¬ 
els to the proper-motion and line-of-sight dispersion profiles together. We compare our dynamical distance 
estimates with literature photometric estimates from the Harris (1996, 2010 edition) globular cluster catalog 
and find that the mean fractional difference between the two types is consistent with zero at just -1.9 ± 1.7%. 

This indicates that there are no signihcant biases in either estimation method and provides an important vali¬ 
dation of the stellar-evolution theory that underlies photometric distance estimates. The analysis also estimates 
dynamical M/L ratios for our clusters; on average, the dynamically-inferred M/L ratios agree with existing 
stellar-population-based M/L ratios that assume a Chabrier initial mass function (IMF) to within -8.8 ±6.4%, 
implying that such an IMF is consistent with our data. Our results are also consistent with a Kroupa IMF, but 
strongly rule out a Salpeter IMF. We detect no correlation between our M/L offsets from literature values and 
our distance offsets from literature values, strongly indicating that our methods are reliable and our results are 
robust. 

Subject headings: globular clusters: individual (NGC104 (47Tuc), NGC288, NGC1851, NGC2808, 
NGC5139 (w Centauri), NGC 5904 (M5), NGC 5927, NGC 6266 (M62), NGC 6341 
(M92), NGC 6388, NGC 6397, NGC 6656 (M22), NGC 6715 (M54), NGC 6752, 

NGC 7078 (M 15)) - proper motions - stars: kinematics and dynamics - stars: distances 
- stars: luminosity function, mass function - globular clusters: general 


1. INTRODUCTION 

Globular clusters (GCs) are among the first objects to have 
formed in the Universe, within 1 Gyr of the Big Bang. They 
are some of the oldest objects for which ages are known and, 
thus, provide crucial constraints on the age of the Universe it¬ 
self (eg. Vandenberg et al. 1996; Carretta et al. 2000; Krauss 
& Chaboyer 2003). As Verde et al. (2013) pointed out, lo¬ 
cal measurements of cosmological parameters (such as the 
age of the Universe) are valuable as they provide cosmology- 
independent tests of favoured cosmological models. 

Three primary methods are used to estimate GC ages: lu¬ 
minosity of the main-sequence turn off (MSTO, eg. Chaboyer 
1995); white-dwarf (WD) cooling curves (Hansen et al. 
2002); and nucleocosmochronology (radioactive age dating of 
stars, eg. Cowan et al. 2002). The MSTO luminosity method 
is the most common, however it requires an accurate determi¬ 
nation of the observed magnitude of the MSTO, accurate red¬ 
denings, and accurate distances. Advances in photometric in¬ 
strumentation have provided both the depth and precision nec¬ 
essary to pinpoint the MSTO within cluster CMDs. Redden¬ 
ing maps, such as those presented in Schlafly & Finkbeiner 
(2011), fulhl the second requirement. However, cluster dis¬ 
tances are still poorly determined and distance uncertainties 
dominate over other sources of error in age determinations. 

Iwatkins @ stsci.edu 
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There are two key points we need to address: 1) we must en¬ 
sure that we are measuring distances accurately; 2) we must 
ensure that different distance measures agree (ie. our methods 
are not biased). 

Cluster distances are also useful to obtain for their own 
merit. They constitute an important rung on the short end of 
the distance ladder, which can again have implications for ver¬ 
ification of cosmological models. Furthermore, distances pro¬ 
vide a third dimension of positional phase space. Many clus¬ 
ters are coincident with substructures in the Galactic halo (eg. 
Bellazzini et al. 2003; Martin et al. 2004; Mackey & Gilmore 
2004); distances to the clusters (and to the substructures) can 
help us determine if there is simply a chance alignment of 
cluster and substructure along the LOS or if they are truly as¬ 
sociated. Evidence of tidal tails has also been found around 
some clusters (eg. Niederste-Ostholt et al. 2010; Kiipper et al. 
2015) and other clusters have been posited as progenitors for 
tidal streams of undetermined origin; distance information 
can help in both cases by providing extra information with 
which to constrain cluster orbits. Cluster orbits, in turn, can 
help to constrain the shape of the inner Galactic halo. 

Cluster distances can be estimated from both photometric 
and kinematic data. Photometric studies typically estimate 
distances via the distance modulus, defined as the difference 
between the apparent and absolute magnitudes of a given ob¬ 
ject or population. Apparent magnitudes are relatively easy 
to obtain. Absolute magnitudes are only well defined for cer¬ 
tain types of objects, often called “standard candles”. Peri¬ 
odic variables, such as Cepheids and RR Lyraes, are common 
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standard candles in the local Universe, although even their 
absolute magnitudes do depend on the metallicity of the star 
and the period of their variability (eg. Catelan et al. 2004). 
Nevertheless, they are very useful distance indicators and are 
widely used; not only for cluster-distance estimation (eg. Del 
Principe et al. 2006), but also for many other science applica¬ 
tions, including to identify substructure in and trace the shape 
of the Galactic halo (eg. Watkins et al. 2009). 

Kinematic studies estimate distances by comparing line-of- 
sight (LOS) velocity dispersions and proper-motion (PM) ve¬ 
locity dispersions. These quantities are measured in different 
units; converting from one set of units to the other requires 
knowledge of the distance. Or, conversely, by assuming that 
the dispersions are equal (or different in some known and 
quantifiable way), the distance can be estimated. Very few 
kinematic distance estimates exist, as only very few clusters 
have measured internal PMs. 

Although, kinematic and photometric distance estimates 
use very different data sets and very different methods, their 
estimates should, of course, be the same for any given clus¬ 
ter. Until now, it has been difficult to assess the agreement 
(or lack thereof) between kinematic and photometric distance 
estimates due to the scarcity of kinematic distance estimates. 
However, in the few cases where both photometric and kine¬ 
matic estimates do exist, they often disagree (eg. Del Principe 
et al. 2006; Watkins et al. 2013). 

Pinpointing the source of the disagreement is complicated 
by the fact that both types of estimate are usually made on a 
cluster-by-cluster basis. So when comparing a set of photo¬ 
metric estimates with a set of kinematic estimates, it is diffi¬ 
cult to determine whether any differences are due to system¬ 
atic biases in either estimation method, or are due to scatter 
from study to study. 

To make progress, we need large, homogeneous sets of 
distance estimates: larger samples will provide more robust 
statistics; homogeneity in the distance estimation process will 
eliminate some potential sources of disagreement so we can 
be confident we are making a meaningful comparison. 

Most GC distances have been measured in individual stud¬ 
ies, however, fortunately, there is a large photometric distance 
analysis that suits our purposes: Harris (1996, 2010 edition, 
hereafter H96) provided a catalogue of cluster properties (in¬ 
cluding distances) for 157 Galactic GCs. No large study has 
yet been done for dynamical distance estimates, as the kine¬ 
matic data has simply not been available. We require both 
LOS and PM data; although both types of data are lacking, 
the dearth of PM data has, so far, been the limiting factor for 
kinematic distance estimates. 

This situation has recently changed. In Bellini et al. (2014, 
hereafter Paper 1), we presented a set of Hubble Space Tele¬ 
scope (HST) PM catalogues for 22 GCs in the Milky Way. 
The catalogues were compiled following a search through 
archival HST data to find fields in Galactic GCs that had been 
previously observed at multiple epochs for different projects. 
The exquisite astrometric sensitivity of HST, coupled with 
baselines of up to 12 years (achievable thanks to the archival 
nature of the project), allowed us to measure PMs with excep¬ 
tional precision. For typical clusters in our sample, the bright, 
well-measured stars have measurement uncertainties of just 
32 /ias yr“* (~L2 km s“' at a distance of 8 kpc, which is typ¬ 
ical for our clusters). Another advantage of these catalogues 
is the sheer volume of data they encompass: the median num¬ 
ber of stars in the full catalogues is ~60,000 and, following 


certain selection cuts as we will describe in Section 2, we use 
a median of 3000 stars per cluster for the analysis we present 
here. 

These catalogues were compiled with two key drivers: 1) to 
facilitate more-detailed studies of each individual cluster; and 
2) for the insight we can gain by considering the population of 
clusters as a whole. There is a wide range of science questions 
that can be addressed both for individual clusters and for the 
whole population, too many to consider in detail here, but see 
the introduction of Paper 1. Here, we seek to advance both 
goals. 

In Watkins et al. (2015, hereafter Paper 2), we presented 
kinematical profiles and maps for the bright stars in each of 
the 22 clusters. Here we build on this analysis by compar¬ 
ing our kinematics with previous LOS studies in order to es¬ 
timate distances for the clusters. This is the first dynamical 
distance study performed for more than a handful of clus¬ 
ters, and greatly increases the number of dynamical distances 
available for Galactic GCs. By comparing these dynamical 
distances estimates with photometric estimates from the liter¬ 
ature, we can assess the consistency of photometric and kine¬ 
matic distance-estimation methods. 

Section 2 introduces the PM catalogues, LOS literature data 
and surface brightness profiles we use for our study. In Sec¬ 
tion 3, we briefly outline our distance estimation method and, 
in Section 4, we use it to estimate distances for our clusters. 
In Section 5, we provide context for our results via compari¬ 
son with previous studies, and look at statistics for the sample 
as a whole. We conclude in Section 6. 

2. CLUSTER DATA 

To estimate cluster distances, we require both PM and LOS 
velocity-dispersion profiles. For our analysis, we use the HST 
PM catalogues described in Paper 1, to which we refer for 
detailed explanations of the data reduction and processing. 
We augment these catalogues with LOS velocity data avail¬ 
able in the literature. Although Paper 1 provides catalogues 
for 22 Galactic GCs, here we are only able to use the clus¬ 
ters for which LOS dispersion profiles have been published 
or for which LOS velocity catalogues have been published 
containing sufficient stars from which to determine a disper¬ 
sion profile. Of the 17 clusters with both PM and LOS data, 
we exclude NGC 6362 and NGC 7099 (M 30) as our PM data 
sets are too small (following the cuts we make below) for our 
distance estimation procedure. 

As we will see in Section 3, we also require a Multi- 
Gaussian-Expansion (MGE)^ fit to the surface brightness pro¬ 
file for each cluster. Fortunately surface brightness profiles 
are available for all of our clusters so we do not need to fur¬ 
ther cut our sample. This leaves us with 15 clusters in total 
for this study. 

Here we briefly describe the PM catalogues, the literature 
LOS data, the surface brightness profiles we use and the MGE 
fitting process. 

2.1. Proper-motion catalogues 

GCs are collisional systems; as the stars undergo a series 
of two-body interactions, they exchange energy and slowly 
move towards a state of energy equipartition. As a result, low- 
mass (faint) stars tend to move faster than high-mass (bright) 
stars. Recently, Trent! & van der Marel (2013) showed that 

^ MGEs are examples of Finite Mixture Models (McLachlan & Peel 2000). 
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clusters never reach full equipartition, but nevertheless, even 
in partial equipartition, we would still expect that stars of dif¬ 
ferent mass will exhibit different kinematic properties. Thus, 
to estimate distances via comparison of PM and LOS disper¬ 
sion profiles, we must ensure that both measurements are for 
stars of similar mass. LOS studies are typically limited only to 
the brightest stars; most data is collected from ground-based 
facilities and only spectra of bright stars achieve sufficient 
signal-to-noise to measure a velocity. As such, we must im¬ 
pose a magnitude cut on our PM sample to select stars of sim¬ 
ilar mass to those observed in LOS studies. 

Furthermore, unbiased kinematics can only be extracted us¬ 
ing stars for which positions (and, hence, PMs) have been 
accurately measured and for which uncertainties have been 
correctly estimated. Poor position measurements - particu¬ 
larly a problem in crowded fields where stars suffer greatly 
from blending - and underestimated error bars will both tend 
to artificially inflate the velocity distributions. The catalogues 
also contain some contaminants, both from the Milky Way 
and from other nearby objects, which must be removed. 

In Paper 2, we described a series of cuts that we make to the 
PM catalogues to ensure that we have a reliable, high-quality 
sample of stars. We refer to Paper 2 for a detailed discussion 
of the cuts, but briefly outline the main points here: 

1. Remove all stars fainter than 1 magnitude below the 
main-sequence turn off. 

2. Remove all stars with Aused/M'ound < 0.8, where Aused is 
the number of data points used for the final calculation 
of the PM and Afound is the total number of data points 
found for each star. 

3. Remove all stars with Tb (DXreduced) > 0.99, where the 
Xreduced values are for the straight-line fits used to de¬ 
termine PMs, Fu is the cumulative distribution function 
for a distribution with D degrees of freedom and, 
here, D = Aused-2. 

4. Remove stars for which the quality of the point-spread 
function (PSF) fit was determined to be low. Quality 
of the PSF fit varies with both magnitude and position 
within the cluster, and also depends on the distance and 
concentration of the cluster. As such, the details of this 
step varied from cluster to cluster. 

5. 3-(T clip in velocity and remove all stars with PM errors 
larger than half the local dispersion. This is an itera¬ 
tive process to ensure that only “good" stars (ie. mea¬ 
surements with enough signal) are used to define the 
boundaries of the final PM and PM error cuts. 

For Paper 2, we required a limited range of stellar masses 
across our sample so that we could neglect the effect of stel¬ 
lar mass on the kinematics. It is this that motivated the initial 
magnitude cut. A cut 1 magnitude below the main-sequence 
turn off was adequate for our purpose of calculating spatial 
kinematic profiles. However, if we are to calculate accurate 
distances via comparison to LOS data, we must be more re¬ 
strictive and impose a brighter magnitude cut. So, we take 
the final samples from Paper 2, extracted via the cuts listed 
above, and make a further magnitude cut, retaining only stars 
brighter than the main-sequence turn off.^ 

^ In principle, we can make this revised magnitude cut at the main- 


2.2. Proper-motion kinematic profiles 

We determine kinematic profiles from our PM data by bin¬ 
ning the cleaned samples in radius and then calculating the 
mean and dispersion of the velocity in each bin. Our PM pro¬ 
files measure the average 1-dimensional PM in the plane of 
the sky, which corresponds to an average over the radial and 
tangential PM components (this is rigorously defined in Sec¬ 
tion 2.6 of van der Marel & Anderson 2010). 

When binning stars, we need to ensure that we have suf¬ 
ficient stars with which to calculate reliable dispersions, but 
we also want to maximise the spatial resolution of our bins. 
As the spatial coverage of the data is highly inhomogenous, 
we And that a variable binning scheme is best suited to our 
needs. Then the kinematics in each bin are determined using 
a maximum-likelihood estimation method that fully accounts 
for biases. Both the binning scheme and the maximum- 
likelihood kinematic estimation are described in full in Paper 
2, to which we refer the reader for further details. 

As described above, the final PM samples we use for this 
paper are different from those we used in Paper 2 as we per¬ 
formed an additional magnitude cut for our present study. As 
a result, the kinematic profiles are not the same as those used 
in Paper 2. 

2.3. Line-of-sight velocity and proper-motion literature data 

Where possible, we use LOS velocity-dispersion profiles 
taken directly from previous studies. However, dispersion 
profiles are not always available; in that case, we use the 
published stellar catalogue of LOS velocities and determine 
our own velocity-dispersion profile. To do this, we proceeded 
similarly as for our PM data: we binned the stars in radius 
(albeit with a simpler binning scheme"^) and used the same 
maximum-likelihood method to determine the kinematics in 
each bin. 

Four clusters also have published PM velocity-dispersion 
profiles; we will later show these profiles in comparison to 
our own but will not use them as part of our distance estima¬ 
tion. All sources of literature data are shown in Table 1 ; for 
each cluster we list the references for the source and the type 
of data taken from the source (LOS dispersion profiles, LOS 
stellar catalogues or PM dispersion profiles). 

We also show the compilations of literature data for each 
of our clusters in Figure 1. Cluster names are shown in the 
bottom-left corner of each panel. Points are coloured accord¬ 
ing to source and data type (LOS or PM), as indicated in the 
legends; colours are different in each panel. Source refer¬ 
ences were abbreviated, the full references are given in Ta¬ 
ble 1. These compilations are also available in Appendix B. 

2.4. MGEfits to surface-brightness profiles 

We begin with the Chebyshev polynomial fits to the V-band 
surface brightness profiles given in Trager et al. (1995). These 
are observed profiles and do not account for interstellar red¬ 
dening. If we only cared about distances, we could proceed 
straight away with the uncorrected profiles, as we only need 

sequence turn off in step 1, however this significantly reduces the number of 
stars in the sample. We then find it difficult to perform the last two stages of 
the cleaning process effectively due to the nature of the cleaning algorithm. 
The two-stage magnitude cut we describe in the text is more effective as it 
limits the mass-range of the initial sample while retaining enough stars with 
which to perform the rest of the cleaning. 

^ For the LOS data sets, we separated the stars into N equally-populated 
bins, where N was chosen separately for each data set; we did not use the 
variable-binning algorithm used for our PM dispersion profiles. 
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Figure 1. Compilations of literature dispersion profiles for each of our clusters. The cluster name is shown in bottom left corner of each panel. The data sources 
and type of data (LOS or PM) are provided in the legend. The sources are abbreviated for convenience; for the full references see Table 1 . PM datasets are shown 
at the distance provided in H96. 
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Table 1 

Sources for literature data. 


Cluster 

(1) 

Source 

(2) 

Type 

(3) 

NGC 104 

Gebhardt et al. (1995) 

LC 


McLaughlin et al. (2006) 

LD, PD 


Lane et al. (2010) 

LD 

NGC 288 

Scarpa et al. (2007) 

LD 


Sollima et al. (2012) 

LD 

NGC 1851 

Scarpa et al. (201 1) 

LD 


Cari'etta et al. (201 1) 

LC 


Liitzgendorf et al. (2013) 

LD 


Lardo et al. (2015) 

LD 


Larson & Seth (priv. comm.) 

LD 

NGC 2808 

Carretta et al. (2006) 

LC 


Liitzgendorf et al. (2012) 

LD 


Lardo et al. (2015) 

LD 


Larson & Seth (priv. comm.) 

LD 

NGC 5139 

van de Ven et al. (2006) 

LD, PD 


Sollima et al. (2009) 

LD 


van der Marel & Anderson (2010) 

PD 


Noyola et al. (2010) 

LD 


Scarpa & Falomo (2010) 

LD 

NGC 5904 

Kimmig et al. (2015) 

LC 

NGC 5927 

Lardo et al. (2015) 

LD 


Simmerer et al. (2013) 

LC 

NGC 6266 

McNamai'a et al. (2012) 

PD 


Liitzgendorf et al. (2013) 

LD 

NGC 6341 

Drukier et al. (2007) 

LC 


Meszaros et al. (2009) 

LC 


Kamann et al. (2014) 

LD 


Kimmig et al. (2015) 

LC 

NGC 6388 

Liitzgendorf et al. (2011) 

LD 


Lanzoni et al. (2013) 

LD 


Lapenna et al. (2015) 

LD 

NGC 6397 

Meylan & Mayor (1991) 

LD 


Gebhardt et al. (1995) 

LC 


Lovisi et al. (2012) 

LC 


Heyletal. (2012) 

PD 

NGC 6656 

Lane et al. (2009) 

LD 


Zloczewski et al. (2012) 

LC 

NGC 6715 

Bellazzini et al. (2008) 

LD 


Ibata et al. (2009) 

LD 


Can'etta et al. (2010) 

LC 

NGC 6752 

Carretta et al. (2007) 

LC 


Sollima et al. (2012) 

LD 


Zloczewski et al. (2012) 

LC 


Lardo et al. (2015) 

LD 

NGC 7078 

Drukier et al. (1998) 

LD 


McNamara et al. (2003) 

PD 


van den Bosch et al. (2006) 

LD,PD 


den Brok et al. (2014) 

LD 


Lardo et al. (2015) 

LD 


Kimmig et al. (2015) 

LD 


Notes. Columns: (1) cluster identification in the NGC catalogue; (2) 
reference for literature data; (3) type of data taken from source (LD = LOS 
velocity-dispersion profile, LC = LOS velocity stellai* catalogue, PD = PM 
velocity-dispersion profile). 


the shape of the luminosity profile for our subsequent anal¬ 
ysis, not the normalisation factor. However, as a by-product 
of our analysis, we will also be able to estimate cluster mass- 
to-light (M/L) ratios and masses; for these M/L and mass es¬ 
timates to be meaningful, we must ensure the surface bright¬ 
ness profiles have been properly corrected. To do this, we use 



Figure 2. MGE fit to the surface-brightness profile for NGC 2808. The black 
line shows the full MGE profile and the coloured lines show the individual 
Gaussian components. 


extinction values Ay from Schlafly & Finkbeiner (2011).^ 
Next, we convert the surface brightness values from 
mag/arcsec^ to L 0 /pc^ via 


S = 




( 1 ) 


where we have used E to denote the surface brightness in units 
of Lq/pc^ and /r to denote the surface brightness in units of 
mag/arcsec^. = 4.83 mag is the V-band absolute mag¬ 
nitude of the Sun and c = 7r/648000 rad arcsec“* is the con¬ 
version factor from arcsec to radians. 

Finally, to perform the MGE fitting, we use the routine de¬ 
scribed in Cappellari (2002).® We start the fitting procedure 
with 15-component fits, however many of the final fits con¬ 
tain fewer Gaussian components as the routine itself discards 
components that make a negligible contribution to the overall 
profile. As an example. Figure 2 shows the surface bright¬ 
ness MGE for NGC 2808, broken down into its constituent 
Gaussians. The coloured lines show the individual Gaussian 
components and the black line shows the total MGE profile. 

The V-band surface brightness profiles and our MGE fits 
are shown in Eigure 3. Both sets of data points are taken from 
Trager et al. (1995): red points show the surface brightness 
data, blue points show their Chebyshev polynomial fits. Our 
MGE fits to the Chebyshev polynomials are shown as black 
lines. 


3. DISTANCE ESTIMATION METHOD 

Dynamical distance estimation relies on the fact that LOS 
velocities are measured as physical distance per unit time, 
while PMs are measured as angular distance per unit time. To 
convert physical distance to angular distance, or vice versa, 
we require the distance to the object. Eor an object with PM ^ 
(in mas yr“'), its transverse velocity v (in km s“') is given by 

v = DCfj,, (2) 

^ In fact, we pulled the extinction values from the NASA/IPAC Extragalac- 
tic Database. 

® While this fitting approach follows the commonly accepted procedure in 
stellar dynamics (e.g. van de Ven et al. 2006; van den Bosch et al. 2008; Em- 
sellem et al. 2011; Liitzgendorf et al. 2012; Feldmeier et al. 2013; Seth et al. 
2014), we note that statistical literature suggests that alternative implemen¬ 
tations (such as maximum likelihood estimation with expectation maximisa¬ 
tion) may be worth exploring in the future. 
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Figure 3. V-band surface brightness profiles for each cluster. The cluster name is shown in the top-right corner of each panel. Red points shows the data compiled 
by Trager et al. (1995). Blue points show the Chebyshev polynomial fits from Trager et al. (1995). Black lines show our MGE fits to the Chebyshev polynomials. 
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where D is the distance of the object (in kpc) and C is a con¬ 
stant that takes care of the remaining unit conversions.^ Dis¬ 
persions (T,, and scale similarly. So, if we knew both cr^ and 
cr,,, we could then infer the distance to the object by finding 
the factor by which the PM velocity dispersion profile must be 
scaled to best match the LOS velocity dispersion. Of course, 
this implicitly assumes that the object is isotropic so that, by 
definition, the PM and LOS dispersions are the same. If the 
object is anisotropic, then the PM dispersion will be some fac¬ 
tor of the LOS dispersion, and this must be accounted for in 
the distance estimation. We showed in Paper 2 that the clus¬ 
ters in our sample are very close to isotropic, particularly at 
their centres, so we will assume that the clusters are isotropic 
over the range for which we have data. 

So by comparing the dispersion profiles we obtain from our 
PM catalogues against LOS velocity dispersion profiles from 
the literature, we can estimate distances to our clusters. How¬ 
ever, our PM dispersion profiles and the LOS dispersion pro¬ 
files are not sampled at the same radii and, even worse, in 
some cases, the radial overlap is small or non-existent. This 
makes it challenging to determine for which distance the two 
profiles best agree. Therefore, instead of comparing profiles 
directly, we fit a simple dynamical model to both dispersion 
profiles and then determine for which distance the best-fit 
models agree. We determine the best-fitting models using 
maximum likelihood analysis. 

Here we describe our dynamical models and then lay out 
the likelihood functions we use to assess the quality of the 
model fits. 


3.1. Models 

We use the Jeans Axisymmetric MGE (JAM) models de¬ 
scribed in Watkins et al. (2013), which were expanded from 
Cappellari (2008). We use only the model-predicted LOS dis¬ 
persion profiles as, under the assumption of isotropy, the LOS 
and PM profiles are the same. The JAM models require pro¬ 
jected luminosity-density and mass-density profiles, both in 
the form of an MGE. In general, these may be different, how¬ 
ever, as we will discuss, we will assume that the mass-density 
profiles have the same shape as the luminosity-density profiles 
so that all that is required is the surface brightness profiles (al¬ 
ready described in Section 2.4) and a mass-to-light ratio. 

These models allow for axisymmetry, anisotropy, rotation 
and a mass-to-light ratio that varies with radius. However, 
for this study, we assume that the clusters are spherical, 
isotropic, non-rotating and have a constant mass-to-light ra¬ 
tio. These assumptions are typical of dynamical cluster stud¬ 
ies (eg. McLaughlin & Meylan 2003; Sollima & Nipoti 2010). 
Although typical, it does not necessarily follow that they are 
correct; however, in general, we believe that these assump¬ 
tions are reasonable for most of the clusters. The excep¬ 
tion is oj Centauri, which we know is highly flattened, shows 
significant rotation and is anisotropic (e.g. van de Ven et al. 
2006; van der Marel & Anderson 2010; Bianchini et al. 2013; 
Watkins et al. 2013) (although we do believe the mass-to-light 
ratio is constant (van de Ven et al. 2006)). Nevertheless, we 
include w Centauri in our analysis, however, we will treat the 
results with caution. We will discuss the model assumptions 
in detail in Section 5.4. 

3.2. Likelihood functions 

^ C = 4.74047 kmyrkpc^* mas“* s“*. 


Let us assume that we have a LOS velocity-dispersion pro¬ 
file consisting of dispersions crL,i±<5L,i (in kms“') measured 
at radii Rl.i (in arcsec) and a PM velocity-dispersion profile 
consisting of dispersions crpj ± (in mas yr“') measured at 
radii Rpj (in arcsec). The units here are important, as we will 
see. 

Let us also assume that we have a dynamical model for 
which we can calculate the dispersion CTmodei (in units of 
kms“') at radius R for a given distance D and mass-to-light 
ratio T^. Under the assumptions we laid out above, both dis¬ 
tance and mass-to-light ratio are the only free parameters in 
our models and both are multiplicative factors: the former be¬ 
cause the models are typically calculated for R in physical 
units whereas we require R in angular units to compare with 
data; and the latter under the assumption that the mass-to-light 
ratio is constant through the cluster. In this case, the model 
dispersion is 


0'model(R,D,T)= y ^ model (R, Do'^o) (3) 

for some reference distance Dq and mass-to-light ratio Tq. 
This simplifies the modelling process as we only need to cal¬ 
culate the model one time for (Do,To) and can then scale the 
model appropriately.® Eor convenience, let us define 


_ ^model (R,Do'^o) 


(4) 


We use likelihood maximisation to determine which models 
best fit the data. Eor a particular distance and mass-to-light 
ratio, the log-likelihood function for fits to the LOS velocity 
dispersion profile is 


ln£L(£>,T) 


E- 


VDTcto (j^Ld) - t7L,i 




(5) 


With LOS velocities alone, the best we can do is determine 
factor 


f=s/DT 


(6) 


that is required to bring the model into good agreement with 
the data. We cannot distinguish how much of / comes from 
the distance scaling and how much comes from the mass-to- 
light scaling. 

Eor the PM velocity dispersion profile, the log-likelihood 
function for a given distance and mass-to-light ratio is 


ln/:p(D,T) = ^-. 


^b{Rpj)ao{Rp.j)-apj] 




(7) 


where the extra factor of DC converts the units of ctq from 
kms“'to masyr“'and b{R) is an anisotropy term that controls 
the size of PM dispersions relative to the LOS dispersions. 
Under the assumption of isotropy, as we use here, 1? = 1 at all 
radii. Again, with PMs alone, the best we can do is determine 
factor 



* As we use V-band surface brightness profiles in our dynamical models, 
all our M/L ratios are in the V-band. 

® In principle, we could use any reference values here; for simplicity, we 
will use reference distance Z)o = 1 kpc and reference mass-to-light ratio T = 
1 Mq/Lq. 
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but cannot distinguish the separate contributions made by dis¬ 
tance and mass-to-light ratio. 

However, with PMs and LOS velocities together, we can 
determine both distance D and mass-to-light ratio T; the total 
log likelihood is then 

In £ (D, T) = In £p (D, T) -I- In £l (£», T) . (9) 

Although D and T are the free parameters of interest, they 
are highly correlated. As / and g are uncorrelated, they are 
more effective to use in a parameter search, as we illustrate in 
Section 4.1. So in practice, we seek to find factors / and g, 
then we determine distance D and mass-to-light ratio T via 

and T = fg. (10) 

g 

4. RESULTS 

Now that we have described both the data and the methods 
we will use, we are ready to estimate distances and M/L ratios 
for our clusters. Using the surface brightness profiles and M/L 
estimates, we can also estimate cluster masses. 

4.1. Distance and mass-to-light ratio estimates 

To recap, we calculate spherical, isotropic, non-rotating 
JAM models for which we assume the mass-to-light ratio 
is constant. Under these assumptions, the model velocity- 
dispersion profiles have only two free parameters - distance 
D and mass-to-light ratio T - and both are multiplicative fac¬ 
tors. We fit these models to the literature LOS and our PM 
velocity dispersion profiles; by combining the two fits we can 
estimate the distance and mass-to-light ratio of the cluster. 

To begin, we perform a simple likelihood maximisation test 
using equation (5) to approximate the best-fitting value of / 
and then a second simple likelihood maximisation test using 
equation (7) to approximate the best-fitting value of g. This 
gives us the region of parameter space where the best model 
is to be found, but does not give an indication of the scatter 
amongst the best-fitting models. To better sample the best-fit 
region, we then use the Markov Chain Monte Carlo (MCMC) 
package EMCEE developed by Foreman-Mackey et al. (2013), 
which is an implementation of the affine-invariant ensemble 
sampler by Goodman & Weare (2010); this approach uses 
multiple trial points (walkers) at each step to efficiently ex¬ 
plore the parameter space. We run our MCMC chain with 
100 walkers for 100 steps, and use the last 20 steps as the 
final sample. 

To determine the distance and mass-to-light ratio for each 
cluster, we first use equation (10) to convert the final MCMC 
sample from (/,g) space to (Z),T) space. Then we take the 
median values of each parameter as the best-fitting values and 
use the 15.9th and 84.1th percentile values as lower and upper 
uncertainty limits as these correspond to the 1 -ct confidence 
interval. 

As an example. Figure 4 shows the final parameter distri¬ 
butions at the end of the MCMC run for NGC 5927. The his¬ 
tograms show the distributions of the individual parameters 

*** These uncertainties reflect scatter in the parameter estimates due to ran¬ 
dom errors; they do not account for any additional sources of scatter that 
could arise if the modelling assumptions were relaxed. Note that we have 
not attempted to constrain the shape of the dispersion profile, we have simply 
used the model profile as a tool to help us match the LOS and PM data and 
estimate distances and M/L values. As such, the quoted uncertainties reflect 
the agreement between the LOS and PM data, not the goodness-of-fit of the 
model to the data. 
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Figure 4. Final parameter distributions at the end of the MCMC run for 
NGC 5927. Histograms show the individual distributions; scatter plots show 
the combined distributions, with points coloured by likelihood from high 
(red) to low (blue). The solid lines show the 50th percentiles the distributions 
and the dotted lines show the 15.9th and 84.1th percentiles of the distribu¬ 
tions. In the upper panels, we show parameters / = y/Wf and g = yjT/D. 
In the lower panels, we show distance D and mass-to-light ratio T directly. 
Comparing the scatter plots, it is clear that D and T are strongly cotTelated 
whereas / and g are uncorrelated, making it more effective to perform our 
parameter search in (/, g) space. 


with a best-fitting Gaussian shown as a black line. The scat¬ 
ter plots show the 2-dimensional distributions of the param¬ 
eters, with the points coloured according to likelihood from 
high (red) to low (blue). The solid lines show the 50th per¬ 
centiles the distributions and the dotted lines show the 15.9th 
and 84.1th percentiles of the distributions. The upper pan¬ 
els show the results for / and g, which show no correlation. 
The lower panels show the results for D and T, highlighting 
their very strong correlation and illustrating why we run our 
MCMC chains in (f,g) space. 

For five of our clusters - NGC 104, NGC 5139, NGC 5927, 
NGC 6388 and NGC 6715 - we restrict the range of the liter¬ 
ature LOS data used for the htting to be only those data points 
that overlap with our PMs. As we discuss in detail in Ap¬ 
pendix A, these cuts are an attempt to mitigate the effect of 
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any inconsistencies outside of the range of our data. 

In Figure 5, we show the results of our distance estima¬ 
tion; the cluster name is shown in the top-right corner of 
each panel. The blue points show the literature LOS disper¬ 
sion profiles used in the fitting process. The red points show 
our PM dispersion profiles, shifted to the best-fitting distance. 
The black lines show model profiles using 100 draws from 
the MCMC sampling. Our distance estimate, with uncertain¬ 
ties, is given in the bottom-left corner of each panel. Where 
available, orange points show literature PM dispersion pro¬ 
files shifted to the best-fitting distance; these are shown for 
comparison only and were not used for the fits. The best val¬ 
ues of /, g and the corresponding distances and mass-to-light 
ratios are given in Table 2; we also give the number of PM 
stars that survived the cleaning and magnitude cuts in Sec¬ 
tion 2 and were used for this analysis. 

The model parameters we estimate here provide the best fit 
to the data under a particular set of assumptions, but this does 
not guarantee that the best-fitting model is actually a good 
fit to the data. A visual inspection of the model fits in Fig¬ 
ure 5 can give an indication of how well the models have per¬ 
formed. There are three main points to consider; the fit of 
the model to the PM dispersion profile, the fit of the model to 
the LOS dispersion profile, and the agreement of the PM and 
LOS dispersions. 

In particular, we are interested in the shape of the model 
dispersion profile compared with the shapes of the data dis¬ 
persion profiles. The only freedom we allow here is to change 
the distance or the (constant) mass-to-light ratio, these param¬ 
eters have no power to change the shape of the model profile. 
Clusters for which the model and data profiles shapes are not 
consistent may indicate clusters for which the assumptions we 
introduced in Section 3.1 break down, as we will discuss fur¬ 
ther in Section 5.4. 

For the majority of the clusters in our sample - NGC 104, 
NGC288, NGC 1851, NGC 2808, NGC5904, NGC6266, 
NGC 6341, NGC 6397, NGC 6656, NGC 6715, NGC 6752 
and NGC 7078 - the agreement between PM dispersions, 
LOS dispersions and model fits appears to be good, that is, 
everything is broadly consistent (albeit with no radial overlap 
between PM dispersions and LOS dispersions for NGC 5904). 
Where available for these clusters, the literature PM data also 
appears to be in good agreement with the rest of the data and 
the models, except for NGC 6397 where the literature PMs sit 
above both LOS data and model at the best-fit distance. 

For the three remaining clusters, the agreements between 
PM data, LOS data and model are less convincing. For 
both NGC 5139 (w Centauri) and NGC 6388, the shape of the 
model profile is a poor fit to both the PM and LOS profiles, 
implying that one or more of our assumptions were not valid 
for these clusters. Nevertheless, despite the poor models, the 
data profiles are consistent with each other so the distance es¬ 
timation seems to have worked successfully and we include 
them in our subsequent analysis. 

For NGC 5927, there is considerable scatter in the best¬ 
fitting models, though the models do seem to trace our PM 
data fairly well. Unfortunately, the restricted radial range of 
the LOS data, leaves us with only 3 data points to use for that 
aspect of the fit. Although not ideal, this is enough for us to 
estimate a distance. Given the good agreement between PM 
data and model here, we also include this cluster in our sub¬ 
sequent analysis. 

Several of the cluster models in Figure 5 show central dips 
in their dispersion profiles. This is not unexpected. Spherical 


isotropic models with steep density profiles generically have 
this feature (eg. Binney & Mamon 1982), and such dips are 
also seen in many dwarf elliptical galaxies (eg. Geha et al. 
2003). Although, if GCs do indeed host IMBHs at their cen¬ 
tres then our models would be underestimates of the disper¬ 
sion near the centre and correctly accounting for the presence 
of an IMBH could make the dips disappear. 

4.2. Mass estimates 

Our primary goal here is to estimate cluster distances. How¬ 
ever, as a by-product of the distance estimation process, we 
estimate the mass-to-light ratio for each cluster. As the sur¬ 
face brightness profiles for the clusters are also known, we can 
use these mass-to-light ratios to get the surface mass-density 
profiles for the clusters and, from there, total mass estimates 
for the clusters. The mass estimates are also shown in Table 2. 

5. DISCUSSION 

We wish to compare our dynamical distance and M/L esti¬ 
mates against previous photometric distances and population- 
synthesis M/L estimates to check for biases and systematic 
offsets. We also consider here the effects of the assumptions 
we made in Section 3.1. 

5.1. Comparison with photometric distance estimates 

We have greatly expanded the number of dynamical dis¬ 
tance estimates available, and all have been estimated using 
the same method. Both of these points will enable us to make 
a meaningful comparison to photometric distance estimates. 

Here, we compare our dynamical-distance estimates with 
the (photometric) estimates from H96, which we include in 
Table 2. These estimates are suitable for this comparison as 
they were all derived in the same way (albeit from different 
datasets). The estimates were made using the mean magni¬ 
tude of the observed horizontal branch Vhb for each cluster. 
A simple calibration relation provides the absolute magnitude 
of the horizontal branch Mv.hb given the cluster metallicity, 
from which the distance modulus p, = Vhb “47v,hb and, hence, 
distance follow. 

In Figure 6, we show the fractional difference between our 
distance estimates and the H96 distance estimates, defined as 
(Doui--DHarris) /Dnarris, as a function of our distance estimates. 
The dotted line marks a difference of zero. The solid line 
marks the mean of the fractional offsets and the shaded areas 
mark the 1-, 2- and 3-tT regions, where a is calculated as the 
error on the mean. The mean and a values are shown in the 
top-right corner of the panel. The most distant cluster in the 
sample, by some margin, is NGC 6715 (M54); the rest of the 
clusters lie within 2 and 12 kpc from the sun. 

The unweighted mean of the offsets is consistent with zero 
at -1.7± 1.9%; this demonstrates that, on average, the dynam¬ 
ical and photometric estimates are in very good agreement 
over the whole sample. Although some clusters do show con¬ 
siderable disagreement between the dynamical and photomet¬ 
ric estimates, this should be further investigated on a cluster- 
by-cluster basis. We see no evidence of a systematic offset 
biasing one method to over- or underpredict compared to the 
other. Furthermore, neither the sign nor the magnitude of the 
offset correlate with cluster distance. 

5.2. Comparison with literature dynamical distance 
estimates 
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Figure 5. Velocity dispersion profiles and model fits for each of our clusters. Our PMs are shown in red, shifted to the best-fitting distance; the literature LOS 
velocities used for the fitting are shown in blue. 100 model draws from the final MCMC sample are shown as black lines. We give the distance estimate with 
uncertainties in the bottom-left comer of each panel; note that the quoted uncertainties reflect the scatter in the fit due to random errors and do not account for 
further systematic errors due to incorrect modelling assumptions. Where available, we show literature PM dispersion profiles (also shifted to the best-fit distance) 
in orange; these are shown for comparison and were not used for the fits. 
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Table 2 

Results: distance and mass-to-light ratio estimates. 


Cluster 

Other names 

N 

/ 

8 

D 

Tv 

^tot 

^literature 

literature 




(kpcM0/L0)O'5 

(M0/L0/kpc)°'^ 

(kpc) 

(M 0 /L 0 ) 

(105 M0) 

(kpc) 

(M 0 /L 0 ) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

NGC 104 

47Tuc 

5631 

2 4i-K).05 
^•^^-0,04 

n ^8 1+0.003 

U.D81-0 004 

4.15 + 0.08 

1.40 + 0.03 

5 57+0-33 
‘^•^'-0.28 

4.5 

2.35 

NGC 288 


713 

4 44+0-24 
-0.21 

0 495+0-011 
U.‘+^J_0012 

9 

^•^■^-0.56 

^•^^-0.10 

0.79tl? 

8.9 

1.97 

NGC 1851 


2810 

3.95 ±0.08 

0.383 + 0.004 

10 32+^'^^ 
i^-7^_0.24 

1.51+0.03 

1.78tl? 

12.1 

1.98 

NGC 2808 


7867 

3 84-^0 05 
-0.06 

0.406 + 0.002 

g 4C+O.I3 

1.56 + 0.02 

5 91+0.22 
‘^•^^-0.25 

9.6 

2.02 

NGC 5139 

(jj Centauri 

16519 

3.72 ±0.05 

0.716 + 0.003 

c iq+0.07 
-^•^^-0,08 

2.66 + 0.04 

34 32+1 "1^ 

7^.7Z-i 43 

5.2 

1.87 

NGC 5904 

M5 

2436 

•1 '50+0.22 

^•^■^-0.25 

0.429 + 0.004 

7 7Q+0.47 
'•^^-0.61 

1 43+0-09 

3.65 + 0.75 

7.5 

1.96 

NGC 5927 


5801 

0 41+0.36 
^•^^-0.38 

0.432 + 0.003 

7 Qi +0.79 
'•^^- 0.88 

1 40+0.15 
l-^°-0.17 

1 44+0.49 
^•^^-0.43 

7.7 

2.93 

NGC 6266 

M62 

7102 

•1 77+0.08 
-0.07 

0.587 + 0.004 

6.42 + 0.14 

2.22 + 0.04 

6 09+0-30 
^•^^-0.33 

6.8 

1.95 

NGC 6341 

M92 

2729 

3.80tl? 

0 426"^-^^ 

U.^ZD_0 005 

0 g3+0.31 

1 62+0-06 

7 8 ^+0-30 
2.57-0 25 

8.3 

1.93 

NGC 6388 


6746 

4.28t!? 

0 392+0-^03 
u.j:?z.-ooo2 

iU.yu_045 

1 6i8+0-06 
^■”°-0.07 

8.27tg 

9.9 

2.55 

NGC 6397 


517 

2.31 + 0.09 

0 967‘^-^22 
U.^O/-0 023 

2 ':5Q+0-13 

7 7'5+0.10 
^•+■^-0.09 


2.3 

1.88 

NGC 6656 

M22 

1924 

/..ji-O 12 

0.816 + 0.010 

2.84 + 0.16 

1 88+0.12 

2 4Q+0.44 
+•^^-0.37 

3.2 

1.87 

NGC 6715 

M54 

3790 

6.61 + 0.10 

0 293‘^-^®- 
u.z,:;j-0 003 

79 C7+0-44 
+^-7/-0,39 

1.94 + 0.03 

11 .83!°:“ 

26.5 

1.88 

NGC 6752 


1860 

2 03+0-06 
^•^■^-0.07 

0.728 + 0.009 

-0.08 

2 14+0.05 
^•^^-0.06 

1.82 + 0.12 

4.0 

1.88 

NGC 7078 

M15 

3107 

3.92 + 0.05 

0.379 + 0.003 

10 36+0' 

1 . 49 + 0.02 

4.95 + 0.19 

10.4 

1.93 


Notes. Columns: (1) cluster identification in the NGC catalogue; (2) alternate name for cluster; (3) number of PM stars used for the distance estimation; (4) / 
parameter used in fitting procedure; (5) g parameter used in fitting procedure; (6) heliocentric distance estimate; (7) V-band mass-to-light ratio estimate; (8) total 
cluster mass estimate; (9) heliocentric distance from H96; (10) V-band mass-to-light ratio from McLaughlin & van der Marel (2005). 
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Figure 6. Comparison of our dynamical distance estimates with literature 
photometric distance estimates from H96. The red points show the fractional 
difference in distance as a function of our distance estimate for each of our 
clusters. The dotted line highlights zero, that is perfect agreement between 
the two estimates. The solid line is the mean of the distance offsets and 
the shaded areas show the 1-, 2- and S-cr regions, where a is the error on 
the mean. The mean and error are also shown in the top-right corner of the 
panel. Overall, the agreement between dynamical and literature distances is 
very good, suggesting that there are no significant biases in either estimation 
method. 


Figure 7. Comparison of dynamical distance and dynamical M/L estimates 
with literature photometric distances (H96) and population-synthesis M/L es¬ 
timates (McLaughlin & van der Marel 2005). The orange points show the 
fractional offsets of our estimates from the literature values; the dotted lines 
highlight zero, or perfect agreement. The unweighted mean and error on the 
mean for the M/L offsets is shown in the top-left comer of the panel. The cor¬ 
responding value of the distance offsets was already show in Figure 6. The 
key point here is that the points are scattered across the whole figure, there is 
no correlation between the direction and magnitude of the distance and M/L 
offsets. 


For five clusters in our sample, dynamical distances have 
previously been estimated, so we now compare our estimates 
to these. As we shall see, overall our estimates are consistent 
with the earlier studies: 

NGC104. McLaughlin et al. (2006) estimated a dynamical 
distance of 4.0 ±0.35 kpc, which is in good agreement with 
our estimate of 4.15 ± 0.08 kpc. 

NGC 5139 (oj Centauri). There have been three dynami¬ 
cal estimates for oj Centauri: 4.8 ±0.3 kpc (van de Ven et al. 
2006); 4.73 ±0.09 kpc (van der Marel & Anderson 2010); 
and 4.59 ±0.08 kpc (Watkins et al. 2013). All three studies 
attempted to account for the inner complexities of the clus¬ 
ter (rotation, anisotropy, flattening), which may explain why 
they are more similar to each other than to our estimate of 
5.19);Qgg kpc (which actually agrees better with the H96 pho¬ 
tometric distance of 5.2 kpc). 

NGC 6266 (M62). McNamara & McKeever (2011) esti¬ 
mated a dynamical distance of 7.054 ±0.583 kpc, which is 
consistent with our estimate of 6.42 ±0.14 kpc at the 1-cr 
level. 

NGC 6397. Heyl et al. (2012) estimated a dynamical dis¬ 
tance of 2.2^g 7 kpc, which is consistent with our estimate of 
2.39!|];Jf kpc.' 

NGC 7078 (M15). There have been two previous dynam¬ 
ical distance estimates for this cluster: 10.3 ±0.4 kpc (van 
den Bosch et al. 2006); and 9.98 ±0.47 kpc (McNamara et al. 
2004). Our estimate of 10.36;|;Q}g kpc is in good agreement 
with both. 

5.3. Comparison with literature mass-to-light ratio 
estimates 

We also wish to compare our mass-to-light ratio estimates 
with those already in the literature. For this comparison, 
we use the mass-to-light ratios from McLaughlin & van 
der Marel (2005), which were estimated from population- 
synthesis models. They adopted cluster metallicities from 
H96 and a common age for all clusters of 13 ± 2 Gyr. Com¬ 


bining stellar population models for a given metallicity and 
age with an assumed initial mass function (IMF), they were 
able to estimate intrinsic colours and average mass-to-light 
ratios. They provided 6 different estimates using different 
stellar-population (SP) codes and assuming different IMFs; to 
begin, we adopt the values calculated using the code presented 
in Bruzual & Chariot (2003) with a Chabrier (2003) disk IMF, 
as they did for their final analysis. These M/L estimates are 
included in Table 2. 

As we did for the distances, we calculate the fractional 
difference between our mass-to-light ratio estimates and the 
McLaughlin & van der Marel (2005) mass-to-light ratio es¬ 
timates, defined as (Tour-TMcLaughiin)/TMcLaughiin- In Fig¬ 
ure 7, we show the fractional-distance differences and the 
fractional mass-to-light differences. The dotted lines mark 
offsets of zero. The unweighted mean M/L offset and error on 
the mean are shown in the top-left corner of the panel. (The 
mean distance offset is already shown in Figure 6.) 

The unweighted mean of the M/L offsets is consistent with 
zero at the 1 .3-a level: -8.8 ± 6.4%. It is clear from the figure 
and from the means and errors that the fractional mass-to-light 
differences are larger, in general, than the fractional distance 
differences. The mass-to-light offsets are typically 20-40%, 
while the distance offsets are typically within 10%. However, 
the scatter in the plot shows no correlation. The direction and 
magnitude of the distance offset is independent of the direc¬ 
tion and magnitude of the M/L offset, thus indicating that our 
method is robust. Moreover, the good agreement of our M/L 
estimates with the McLaughlin & van der Marel (2005) esti¬ 
mates implies that a Chabrier (2003) disk IMF is consistent 
with our data for all clusters. 

5.3.1. Initial mass functions 

As previously mentioned, McLaughlin & van der Marel 
(2005) provided 6 different M/L estimates using different SP 
codes and different IMFs. We now compare our dynamical 
M/L estimates with each of these population-synthesis esti- 
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Table 3 

Consistency of different population-synthesis M/L estimates. 


Model 

(1) 

Code 

(2) 

IMF 

(3) 

AT 

(4) 

1 

BC 

Chabrier (disk) 

-8.8 ± 6.4 

2 

BC 

Salpeter 

-45.7 ± 3.9 

3 

PEGASE 

Chabrier (disk) 

-17.9 ± 5.8 

4 

PEGASE 

Salpeter 

-50.6 ± 3.6 

5 

PEGASE 

Chabrier (GC) 

-8.6 ± 6.6 

6 

PEGASE 

Kroupa 

-15.0 ± 6.1 


Notes. Columns: (1) population-synthesis model; (2) SP code used to 
compute model (BC=Bruzual & Chariot 2003, PEGASE=Fioc & 
Rocca-Volmerange 1997 v2.0); (3) IMF used for model (sources: Chabrier 
2003; Salpeter 1955; Kroupa et al. 1993); (4) fractional offset of the 
population-synthesis M/Ls from our dynamical M/Ls. The 
population-synthesis M/L estimates were all taken from McLaughlin & van 
der Marel (2005). 

mates. Once again, we calculate the fractional offset of our 
measurements from the population-synthesis estimates and 
calculate the mean and the standard error on the mean. The 
models are summarised and M/L offsets are reported in Ta¬ 
ble 3. 

Irrespective of the SP code used to perform the modelling, 
we see that our dynamical M/L estimates and the population- 
synthesis estimates using a Salpeter (1955) IMF (models 2 
and 4) are significantly different. The offsets are at the level 
of 12-a using the Bruzual & Chariot (2003) code and 14-cr 
using the PEGASE code (Fioc & Rocca-Volmerange 1997). 
We strongly rule out a Salpeter (1955) IMF for Galactic GCs. 

The remaining models used a Chabrier (2003) disk IMF 
(models 1 and 3), a Chabrier (2003) GC IMF (model 5) or 
a Kroupa et al. (1993) IMF (model 6). All are consistent with 
our M/L estimates, given the random errors inherent in our 
analysis and uncertainties due to the modelling and details of 
the SP codes. 

5.3.2. Mass-to-light ratio variation with metallicity 

Strader et al. (2011) estimated dynamical V-band M/L ra¬ 
tios for a set of 200 globular clusters in M31 and found that 
their M/L estimates tend to decrease with increasing cluster 
metallicity. This behaviour directly opposes predictions from 
population-synthesis models that M/L ratios should increase 
with increasing cluster metallicity. So let us now consider 
how our M/L estimates vary with cluster metallicity, which 
values we take from H96. 

In Figure 8, we show our dynamical M/L estimates as a 
function of cluster metallicity (blue points). For comparison, 
we show the Strader et al. (2011) M31 cluster data (black 
points) and the McLaughlin & van der Marel (2005) M/L esti¬ 
mates for our clusters (green points, again using metallicities 
from H96). The distribution of our M/L estimates is consis¬ 
tent with the Strader et al. (2011) estimates across the full 
range of metallicities and exhibits the same decrease in M/L 
that Strader et al. (201 1) saw at the metal-rich end of the M31 
cluster distribution. Our estimates are also consistent with the 
McLaughlin & van der Marel (2005) population-synthesis es¬ 
timates for the metal-poor clusters (as we would expect given 
the generally good agreement overall discussed in the previ¬ 
ous section) but deviate for the most metal-rich clusters. 

Strader et al. (201 1) argued that this behaviour is not caused 
by dynamical evolution inside the cluster and instead favour 
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Figure 8. M/L estimates as a function of metallicity. The blue points show 
our dynamical M/L estimates, with metallicities taken from H96. For com¬ 
parison, the black points show the dynamical M/L estimates for a sample of 
M31 GCs from Strader et al. (2011), which tend to decrease with increasing 
metallicity; the green points show the population-synthesis M/L estimates 
from McLaughlin & van der Marel (2005) for our clusters, which tend to 
increase with increasing metallicity. Our estimates are consistent with the 
Strader et al. (2011) M31 clusters over the full range of metallicities but de¬ 
viate from the population-synthesis estimates at high metallicities. 

a modification to population-synthesis models (in the form 
of a shallower mass function for metal-rich GCs) to bring 
the dynamical estimates and population-synthesis estimates 
into better agreement. However, Shanahan & Gieles (2015) 
showed that this phenomenon may indeed be dynamical in 
origin and attribute it to mass segregation in the cluster. It is 
beyond the scope of this work to investigate further the ori¬ 
gin of this phenomenon, however the full PM catalogues are 
well-suited to studying mass segregation in clusters and this 
is a topic we plan to address in future papers. 

5.4. Effect of modelling assumptions 

For the modelling part of this analysis, we assumed that 
our clusters are spherical, isotropic, non-rotating and have a 
constant mass-to-light ratio. For 12 of our 15 clusters, the 
models were generally good fits to both the PM and LOS data. 
This suggests that our assumptions were reasonable for these 
clusters; we cannot be certain, though it seems convenient that 
the effect of a number of incorrect assumptions would cancel 
out so well. However, let us briefly consider what we do know 
about these clusters and the quality of the assumptions we 
have made. 

5.4.1. Shape 

The assumption of sphericity implies that quantities of in¬ 
terest vary only with radius and not with position angle so we 
can consider 1-dimensional radial profiles. We demonstrated 
in Paper 2 that many of our catalogues contain sufficient data 
to study 2-dimensional spatial variations in velocity disper¬ 
sion, however most LOS velocity datasets are too small to 
obtain 2-dimensional dispersion maps. (Indeed, many LOS 
velocity datasets are so small that obtaining 1-dimensional 
profiles can be a challenge.) Furthermore, we require sur¬ 
face brightness profiles for our modelling here, and these are 
only available as 1-dimensional profiles for most clusters. De¬ 
riving 2-dimensional surface-brightness profiles is beyond the 
scope of this paper. Reassuringly, the 2-dimensional disper- 
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sion maps we showed in Paper 2 did not reveal a high degree 
of flattening in our clusters. 

White & Shawl (1987) provide ellipticities for 14 of the 
well-fit clusters (only NGC 288 is missing from their list): the 
values range from 0.01 for NGC 6266 to 0.14 for NGC 5904 
and NGC 6656. In fact, there are four clusters with elliptic¬ 
ities of 0.1 or higher. We might expect the assumption of 
sphericity to break down for these moderately flattened clus¬ 
ters, however we see no correlation between distance offset 
and ellipticity. This supports the notion that ignoring elliptic- 
ity does not systematically bias the distance estimates. 

5.4.2. Anisotropy 

We showed in Paper 2 that the clusters in our sample are 
very close to isotropic at their centres. Some clusters do show 
some radial anisotropy with increasing radius, but the effect 
is very mild. In most cases, the LOS velocity data extends 
beyond our PM data, however, even if we extrapolate the 
anisotropy profiles from Paper 2 to cover the range of the LOS 
velocity data, the anisotropy is still small. Paper 2 measured 
projected anisotropy and not intrinsic anisotropy, which is key 
when considering LOS data alongside our PM data. However, 
Section 2.7.3 of van der Marel & Anderson (2010) shows that 
projected and intrinsic anisotropy are directly linked, so our 
expectation of isotropy or, at most, mild anisotropy still holds. 

Furthermore, we average out the impact of anisotropy by 
combining both radial and tangential PMs together when cal¬ 
culating our PM dispersion profiles. For steep density pro¬ 
files, anisotropy will tend to increase one component of mo¬ 
tion while decreasing another component of motion by an ap¬ 
proximately equal factor so that their mean remains approxi¬ 
mately the same. Although this behaviour can break down for 
flatter density profiles, such as we observe for GCs (eg. Cap- 
pellari 2015); we performed tests to verify that this is not the 
case for our clusters. In general, averaging the PMs together 
does tend to wash out any anisotropy. 

Finally, the good agreement between our distances and the 
Harris photometric estimates suggests that the assumption of 
isotropy is reasonable. 

5.4.3. Rotation 

We assume that the clusters are not rotating - this is a nec¬ 
essary assumption because our PM catalogues contain only 
relative motions and not absolute motions, thus washing out 
any rotation signatures. As we discussed in detail in Paper 2, 
the lack of any rotation signature does not affect the disper¬ 
sions that we measure. We also know that rotation will gen¬ 
erally increase with radius and, as our PM measurements are 
confined to the central regions, we do not expect significant 
signal from rotation in our PM data. 

For the LOS data sets, rotation has the potential to be more 
of a concern. The key point to consider here, again, is that we 
expect little to no rotation at the centre where the LOS data 
overlaps with our PM data. The good agreement between 
LOS data and model in the outer regions for most clusters 
would suggest that rotation is not significant in most cases. 
However, recall that there were 5 clusters for which we re¬ 
stricted the range of the LOS data to be only the region that 
overlapped with our PM data due to poor model fits. For these 
clusters, it is possible that rotation may become significant 
in the outer regions and should be included for an improved 
model fit outside of the range of our PM data. 

5.4.4. Mass-to-light ratio 


We assume that mass follows light in clusters, so the mass- 
density profile is simply a scaled version of the luminosity- 
density profile; the scale factor is, of course, the M/L ratio. As 
they evolve, clusters move towards a state of energy equipar- 
tition, the effect of which is that massive (bright) stars tend to 
move more slowly than the low-mass (faint) stars. As a re¬ 
sult, the more-massive (brighter) stars sink towards the centre 
of the cluster and, therefore, have a more concentrated den¬ 
sity profile than the less-massive (fainter) stars. This is mass 
segregation. Owing to the different luminosity-density and 
mass-density profiles created by this mass segregation (Krui- 
jssen 2008), we might expect that the cluster M/L will change 
with radius. 

In principle, our dynamical models could allow the M/L 
to change as a function of radius, however, for our analysis, 
we assume that the mass-to-light ratio is constant through the 
cluster. The M/L at a fixed radius will affect both PM and 
LOS measurements equally, so a radially-varying M/L is only 
an issue in cases for which the radial coverage of the PM data 
and LOS data are different. Once again, the good agreement 
between model and data, even in cases where the radial cover¬ 
age of the different data types are not consistent, indicates that 
the assumption of constant (or nearly constant) M/L through 
the clusters is reasonable. It is possible that the clusters for 
which we restricted the range of the LOS data to have an M/L 
that varies with radius, but the effect is mitigated by consider¬ 
ing only the LOS data that overlaps with our PM data. 

Further support for the assumption of a constant M/L comes 
from van de Ven et al. (2006) who showed that the mass-to- 
light ratio of oj Centauri (NGC 5139) is constant over a wide 
range of radii. Though it is worth noting that ui Centauri may 
harbour an intermediate-mass black hole (IMBH) at its centre 
- results are conflicting (van der Marel & Anderson 2010; 
Noyola et al. 2010). IMBHs are believed to limit the amount 
of mass segregation in a cluster because they act as a source of 
energy at the cluster centre (Baumgardt et al. 2005; Gill et al. 
2008); so if w Centauri does host an IMBH, we would expect 
little mass segregation and a fairly flat mass-to-light ratio. 

IMBHs are not the only mechanism by which mass segre¬ 
gation can be quenched; any source of energy at the centre 
will achieve the same effect. A significant population of bi¬ 
naries or a segregated population of stellar-mass black holes 
at the centre of the cluster can also inhibit mass segregation 
throughout the rest of the cluster (eg. Mackey et al. 2008). 
So there are a number of factors that determine the degree of 
mass segregation and, hence M/L variation in a cluster. 

The degree of segregation within clusters has been studied 
for individual clusters (eg. Pasquato et al. 2009; Beccari et al. 
2010) and for large samples of clusters (eg. Goldsbury et al. 
2013). 14 of our 15 clusters are in the Goldsbury et al. (2013) 
sample; most show some mass segregation but it is relatively 
mild. A notable exception is NGC 6388, which shows consid¬ 
erable mass segregation. This may explain why the shape of 
the Jeans models deviated from the shape of the velocity dis¬ 
persion profiles, though, recall that the velocity profiles them¬ 
selves did show good agreement, despite the poor shape of the 
model. 

6. CONCLUSIONS 

We have determined dynamical distance and M/L estimates 
for 15 Galactic globular clusters. To do this, we combined 
PM dispersion profiles determined from the HST PM catalogs 
described in Paper 1 with LOS dispersion profiles compiled 
from literature sources. To these dispersion profiles, we then 
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fitted dynamical models from which we were able to estimate 
distances and M/Ls. For the majority of our clusters, no PM 
data was previously available and so this is the first dynamical 
distance and M/L study performed. 

Although the distance estimates themselves are a useful 
product of this work, the driving force behind this study was 
to check for consistency between dynamical distance esti¬ 
mates and photometric distance estimates. To this end, in 
Section 5.1, we compared our dynamical distances with pho¬ 
tometric distance estimates from H96. A key point in this 
comparison is the homogeneity of the data analyses that we 
compared; the data for each cluster is necessarily different as 
each cluster is unique, however, both here and in H96, the 
data was treated in a uniform way. This helps to eliminate 
some possible sources of contention and allows for a more 
robust comparison. 

We showed that the dynamical and photometric distance es¬ 
timates were consistent with each other within their uncertain¬ 
ties. The mean fractional difference between the two types of 
estimate was just -1.9 ± 1.7%. This indicates that there are 
no significant biases in either method that may cause one to 
under- or over-estimate the cluster distances with respect to 
the other. This also provides an important validation of the 
stellar evolution theory that underlies photometric distance es¬ 
timates. Furthermore, neither the sign nor the magnitude of 
the fractional distance offset correlate with cluster distance. 

A side product of this analysis was the estimation of dy¬ 
namical M/L values for each cluster. We compared these 
against the population-synthesis M/L values determined in 
McLaughlin & van der Marel (2005). Pleasingly, the val¬ 
ues agreed within 1.3 -ct (the mean fractional M/L offset was 
-8.8 ±6.4%), and we detected no correlation between our 
M/L offsets from literature values and our distance offsets 
from literature values. This provides a strong indication that 
our methods were reliable and our results are robust. 

Further, we compared our dynamical M/L estimates with 
a set of population-synthesis M/L estimates (also from 
McLaughlin & van der Marel 2005) which assumed differ¬ 
ent IMFs: we ruled out a Salpeter (1955) IMF at high sig¬ 
nificance, but were unable to distinguish between a Chabrier 
(2003) disk IMF, a Chabrier (2003) GC IMF and a Kroupa 
et al. (1993) IMF. We also found that our M/L estimates tend 
to decrease with increasing metallicity; this behaviour directly 
opposes predicts from SP models, but supports a similar trend 
observed in a sample of 200 M31 clusters by Strader et al. 
( 2011 ). 

Paper 1 provided catalogues for 22 Galactic globular clus¬ 
ters, however we were only able to use 15 clusters for this 
analysis. Of the clusters we did not include here, two clus¬ 
ters - NGC 6362 and NGC 7099 - have LOS literature data, 
however we found that our PM catalogs (after cleaning) were 
too small to perform a reliable analysis. For these, it is possi¬ 
ble that a more tailored cleaning algorithm would yield more 
usable stars than the one-size-fits-all cleaning we performed 
here - this is beyond the scope of this paper but would be an 
interesting topic for future work. 

For the remaining 5 clusters - NGC 362, NGC 6441, 
NGC 6535, NGC 6624 and NGC 6681 - there was no LOS 
data available in the literature that was sufficient for our pur¬ 
poses here. This illuminates an important shift in dynamical 
cluster studies. Although LOS measurements and PM mea¬ 
surements are both lacking for clusters, until now, lack of PM 
data has been biggest drawback. However, we are now in a 
situation where we have some clusters with LOS data but no 


PM data, and some clusters with PM data but no LOS data. 
These five clusters would be prime targets for observing cam¬ 
paigns in the near future in order to obtain larger LOS data 
sets that could augment the kinematic information provided 
by our PM catalogues. 
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APPENDIX 

RESTRICTING LITERATURE DATA 

During preliminary analysis of our clusters, we found that 
five had several inconsistencies between our PM data, liter¬ 
ature LOS data, and the model fit. Eigure 9 shows the ini¬ 
tial results for these five clusters. As for Eigure 5, red points 
show our PM data and blue points show literature LOS data; 
the black line show the models fits to the red and blue points. 
Orange points show literature PM data, where available, for 
comparison only and were not used for fitting. Let us con¬ 
sider each in turn. 

NGC 104. Here, the model is a reasonable fit to the LOS 
data where it overlaps with our PM data, but underestimates 
the LOS data in the outer regions. This is likely due to the 
break down of our assumptions in the outer regions of the 
cluster (see discussion in Section 5.4). The anomalous two 
outermost data-points are attributable to evaporation (Lane 
et al. 2010); their inclusion or exclusion does not significantly 
alter the fit. 

NGC5139. Here, the models do a very poor job at fitting 
both the LOS dispersions and the PM dispersions. However, 
as we mentioned in Section 3.1, w Centauri is well-known 
to be anisotropic, rotating and highly-flattened, so it is not 
surprising that even the best model does not fit the data well. 
We expect that anisotropy, rotation and flattening will all be 
more significant and intermediate and large radii, that is, in 
regions beyond the extent of our PM data. 

NGC5927. The peculiar shape of the surface brightness 
profile for this cluster (see Eigure 3) gives rise to a peculiar 
shape for the model dispersion profiles, though the models do 
seem to trace our PM data fairly well. However, the models 
fail to accurately reproduce the shape of the LOS dispersion 
profile outside of the range of our PM data. 

NGC 6388. Here, the model is a very poor fit to the data 
in the central regions of the cluster. This is due to the con¬ 
siderable controversy regarding the shape of the dispersion 
profile near the centre. Studies of individual radial velocities 
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Figure 9. Full fits for restricted clusters. 


(Lanzoni et al. 2013; Lardo et al. 2015) argue that the disper¬ 
sion prohle stays flat at the centre, while studies using Inte¬ 
gral Field Spectroscopy (IFS) measure a dispersion that rises 
steeply towards the centre Lutzgendorf et al. (201 1) (as would 
be expected if the cluster harbours an intermediate-mass black 
hole at its centre). The PM profile used for the present study 
does not extend far enough into the centre to comment on ei¬ 
ther of these claims.*^ 

NGC 6715. For this cluster, the model fails to trace the cen¬ 
tral rise of the literature LOS dispersion profile or the sharp 
upturn in the outer regions of the cluster. However, the model 
is broadly consistent with the LOS data in the region of over¬ 
lap with our PM data. 

In an attempt to improve the fits, we restrict the range of the 
literature LOS data used for the fitting of these clusters to be 
only those data points that overlap with our PMs. That is, if 
our PM dispersion profile has an innermost bin at ± 
and an outermost bin at /?out ± A/?outi we restrict the range of 
the literature data to be — Af?in < R < Af^out- 

For some of these clusters, the inconsistencies are only ap¬ 
parent outside of the range of our PM data, so by restricting 
the data in this way, we can eliminate problem areas. For oth¬ 
ers, the fits seem generally poor, but by restricting the range 
of the LOS data to the region of overlap with our PM data, 
we hope to at least mitigate any problems outside of the re¬ 
gion of our data. As discussed further in Section 4.1, the 
fits are greatly improved for NGC 104 and NGC 6715, but 
still show some inconsistencies for NGC 5139, NGC 5927 and 
NGC 6388. 

It is important to note that restricting the LOS data as we 
have described here is not the same as fudging the result. In 
essence, our goal is simply to align the LOS and PM profiles 
in the radial range where they overlap. Employing dynamical 
models, as we have done here, allows us to also use the data 
in the radial range where they do not overlap. This is only 

This will require more careful cleaning of the data and more detailed 
modelling, and will be the topic of a future paper. 


improvement if the models are to be trusted and provide a 
good fit. If they do not fit the data at large or small radii, then 
this approach actually makes the results worse than when only 
using the data in the radial overlap region. 


COMPILATION OF LITERATURE DATA 

Table 4 provides the compilation of literature velocity dis¬ 
persion profiles used for this paper, as described in Section 2.3 
and summarised in Table 1 . This table contains profiles for all 
clusters together. The columns list; the name of the cluster; 
the radial coordinate (and uncertainty, where available); the 
velocity dispersion, with uncertainty, of the stars in the bin; 
the type of data (PM or LOS); and the data source. 


REFERENCES 


Astropy Collaboration, Robitaille, T. R, Tollerud, E. J., et al. 2013, A&A, 
558, A33 

Baumgardt, FI., Makino, J., & Hut, P. 2005, ApJ, 620, 238 
Beccari, G., Pasquato, M., De March!, G., et al. 2010, ApJ, 713, 194 
Bellazzini, M., Ibata, R., Ferraro, F. R., & Testa, V. 2003, A&A, 405, 577 
Bellazzini, M., Ibata, R. A., Chapman, S. C., et al. 2008, AJ, 136, 1147 
Bellini, A., Anderson, J., van der Marel, R. P, et al. 2014, ApJ, 797, 115 
Bianchini, P., Varri, A. L., Berlin, G., & Zocchi, A. 2013, ApJ, 772, 67 
Binney, J., & Mamon, G. A. 1982, MNRAS, 200, 361 
Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 
Cappellaii, M. 2002, MNRAS, 333, 400 
—. 2008, MNRAS, 390,71 
—.2015, ArXiv e-prints 

Carretta, E., Bragaglia, A., Gratton, R. G., et al. 2006, A&A, 450, 523 
Carretta, E., Bragaglia, A., Gratton, R. G., Lucatello, S., & Momany, Y. 
2007, A&A, 464, 927 

Carretta, E., Gratton, R. G., Clementini, G., & Fusi Pecci, F. 2000, ApJ, 533, 
215 

Carretta, E., Lucatello, S., Gratton, R. G., Bragaglia, A., & D’Orazi, V. 

2011, A&A, 533, A69 

Carretta, E., Bragaglia, A., Gratton, R. G., et al. 2010, A&A, 520, A95 
Catelan, M., Pritzl, B. J., & Smith, H. A. 2004, ApJS, 154, 633 
Chaboyer, B. 1995, ApJ, 444, L9 
Chabrier, G. 2003, PASP, 115, 763 

Cowan, J. J., Sneden, C., Buries, S., et al. 2002, ApJ, 572, 861 
Del Principe, M., Piersimoni, A. M., Storm, J., et al. 2006, ApJ, 652, 362 
den Brok, M., van de Ven, G., van den Bosch, R., & Watkins, L. 2014, 
MNRAS, 438, 487 

Drukier, G. A., Cohn, H. N., Lugger, P. M., et al. 2007, AJ, 133, 1041 











HSTPROMO GCs. III. Dynamical distances and mass-to-light ratios 


17 


Table 4 

Compilation of literature dispersion profiles. 


Cluster 

R 

AR 

cr 

Acr+ 

Act 

type 

source 


(arcsec) 


(km s“^) 




(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

NGC 104 

13.525 

4.593 

11.967 

0.939 

0.939 

LOS 

Gebhardt et al. (1995) 


24.967 

2.510 

12.898 

1.079 

1.079 

LOS 

Gebhardt et al. (1995) 


32.605 

2.217 

12.582 

1.096 

1.096 

LOS 

Gebhardt et al. (1995) 


40.639 

2.227 

12.564 

0.839 

0.839 

LOS 

Gebhardt et al. (1995) 


48.045 

2.293 

10.726 

0.842 

0.842 

LOS 

Gebhardt et al. (1995) 


57.730 

A.lll 

10.309 

0.796 

0.796 

LOS 

Gebhardt et al. (1995) 


140.017 

47.826 

10.512 

0.744 

0.744 

LOS 

Gebhardt et al. (1995) 


457.745 

157.683 

8.015 

0.643 

0.643 

LOS 

Gebhardt et al. (1995) 


7.500 

0.000 

9.760 

1.640 

2.090 

LOS 

McLaughlin et al. (2006) 


15.000 

0.000 

11.820 

1.050 

1.130 

LOS 

McLaughlin et al. (2006) 


Notes. Columns: (1) cluster ID; (2-3) radius and uncertainty; (4-6) dispersion and uncertainty; (7) data type; (8) data source. 
(This table is available in its entirety in a machine-readable form.) 
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